Survey of Quadratic Programming Solvers for Applying the Unilateral Incompressibility Constraint in Pressure Solving

In a continuum approach to the simulation of granular materials self intersection between particles becomes an issue. Even with a solution to incompressible flow, the incompressibility constraint won’t prevent particles from intersecting each other. The unilateral incompressibility constraint [1] provides a constraint to prevent most of these self-intersections that are visually obvious and can significantly reduce visual efficacy of a simulation. This article will explore the performance of a variety of solvers to solve the linear problems generated by my choice of discretization.